Scaling properties of the Anderson model in the Kondo regime studied by aGaW 

formalism 
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The symmetric Anderson model for a single impurity coupled to two leads is studied at strong 
interaction using the GW approximation within the aGaW formalism. We find that the low energy 
properties show universal scaling behavior in the asymptotic regime. While the GW scaling functions 
are similar in form to the scaling functions known from the numerically exact solution, they are 
characterized by a different parameter value indicating that GW fails to describe correctly spin 
correlations between the impurity and lead electrons. We also compare the GW and exact Kondo 
scales for a broad range of the interaction strength. In contrast to the exponential behavior shown 
by the exact solution, the GW Kondo scale depends algebraically on the interaction strength. 



I. INTRODUCTION 

As the size of modern electronic components is pushed 
toward the molecular limit, the effects of electron- 
electron interactions can be dramatically enhanced due 
to quantum confinement. Therefore, understanding 
strongly correlated phenomena in nanoscale electron 
transport is important for the development of future elec- 
tronic devices. In this regard, the Anderson impurity 
model 1 has been very useful, being able to capture the 
esential physics of complex transport phenomena. Exam- 
ples include magnetic nanojunctions composed of atoms, 
molecules or quantum dots 2 - - — that show a conductance 
minimum in the linear-response at low temperature, simi- 
lar in origin to the well-known resistance minimum shown 
by metals with magnetic impurities: as temperature is 
lowered, the conductance decreases due to suppressed 
charge fluctuations (the Coulomb blockade effect), and 
then increases due to enhanced spin-flip scattering (the 
Kondo effect). 

Presently these strongly correlated phenomena can- 
not be described entirely from first principles. Recent 
theoretical approaches^ rely on ab initio, perturba- 
tive methods to extract parameters for a specifically 
designed Anderson model that is then solved by non- 
perturbative means such as the numerical renormaliza- 
tion group (NRG) method. Using entirely ab initio ap- 
proaches to gain fundamental understanding of strongly 
correlated electron trasport remains an open challenge 
even for idealized nanojunctions. 

The spin-1/2 single-impurity Anderson model (SIAM) 
is a basic model in condensed-matter physics and a 
paradigm for a system of strongly interacting electrons. 
At equlibrium, its thermodynamic properties can be ob- 
tained exactly by means of the Bcthe-Ansatz, while its 
dynamic (excited state) properties can be obtained nu- 
merically with desired precision, using a variety of non- 
perturbative methods such as NRG and quantum Monte 
Carlo (QMC). At non-equilibrium, a major step toward 
solving the symmetric SIAM has been very recently 
achieved via diagrammatic QMCS. 

Considering the SIAM as an idealized nanojunction, it 
is tempting to study its properties in the strong inter- 



action regime using a many-body perturbative approach 
that is amenable to ab initio calculations. Progress along 
this direction could open the door for complete first- 
principles studies of strongly correlated phenomena in 
realistic nanojunctions. That such progress is possible is 
indicated by an important property of the SIAM, namely 
that the system is always in the Fermi-liquid state as it 
does not undergo any phase transition; thus, perturba- 
tion expansion about the impurity Coulomb interaction 
U is guaranteed to converg o 10 ' 11 . 

Treating the dynamic properties of the SIAM with low- 
order (up to fourthi^r— ) perturbation expansion about U 
has provided so far only a qualitative description of the 
interesting features that the impurity density of states 
develops with increasing U (such as the emergence of the 
sharp Kondo peak near the Fermi level Ep or the forma- 
tion of Hubbard bands away from Ep near ±17/2), and 
then only for not too large interaction strength. In par- 
ticular, low-order expansions in U yield a non-universal 15 
Kondo resonance which narrows algebraically 16 in U in- 
stead of exponentially. 

An alternative approach that we consider involves per- 
turbation expansion about the screened Coulomb interac- 
tion W— . The first order perturbation expansion about 
W, known as the GW approximation 15 (GWA), is one 
of the best available tools in describing the quasipar- 
ticle properties of real materials^— characterized by 
not too strong electron correlation. The GWA has been 
applied also to the SIAM-- 4 "? 7 -, showing clearly that it 
fails to reproduce the Hubbard bands unless it artificially 
breaks spin symmetry. However, an important question 
rcmain a 25 ' 26 : how well does GW describe the Kondo 
peak? That this issue is still open is due in part to the 
fact that a non-magnetic GW solution has been elusive 
in the strong interaction regime near half-filling and low 
temperature and bias voltage 2 ^. Our work answers the 
above question by calculating such GW solution in the 
electron-hole symmetric case, where charge fluctuations 
are most suppressed. 

An important conclusion of our work is that GW shows 
universal scaling behavior in the Kondo regime. A care- 
ful analysis of this behavior shows that GW neither de- 
scribes satisfactorily the Kondo resonance nor the ther- 
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modynamic and transport properties that depend on it. 
On a more positive note, we have identified, employing 
the recently developed aGaW formalism 2 ^ a new GW 
flavor which shows marked improvement over previously 
used GW flavors; this finding might be useful in applica- 
tions of GW to moderately correlated electron phenom- 
ena. 

The rest of the paper is organized as follows. The 
model Hamiltonian and the underlying aGaW formal- 
ism are described in Sec. II. Section III presents results 
of GW calculations for the in- and out-of-equilibrium low- 
energy properties of the SIAM in the strong interaction 
regime, compared (at equilibrium) to results known from 
the numerically exact or Bethe-Ansatz solution. Simple 
arguments attempting to relate the parameters that con- 
trol the GW scaling functions and the Kondo scale are 
provided in the Appendix. 

II. FORMALISM 

A. Hamiltonian 

We consider the Anderson modei^ for an impurity cou- 
pled symmetrically to (left and right) non-interacting 
semi-infinite leads. The non-interacting part of the 
Hamiltonian describing the spin- 1/2 impurity has the 
form: 

H imp = ^2 e dd\d a (1) 

where d\ creates an electron with spin a on the impurity. 

We describe the leads within the infinite band-width 
limit (this approximation has negligible effect on low- 
energy quantities that are the subject of our paper) . The 
hybridization function describing the coupling between 
leads and impurity is then just a constant, the effective 
coupling strength that we denote by A. 

The electron-electron interacting part of the Hamilto- 
nian is localized on the impurity and represented by the 
usual U term: 

H c . c = Un d1: n di = i ^ rfa d p V <xa',PP'd0'da' ( 2 ) 

where rida is the electron occupation number of spin a on 
the impurity and U is the repulsive Coulomb interaction 
between spin up and spin down impurity electrons. 

There are several physically appealing choices for the 
two-particle electron-electron interaction V that give the 
same H e . e in Eq. @. In this work, we consider two dif- 
ferent spin-dependent (and free of self-interaction effects) 
versions: 

Vl«,w = -u(i - M(! - <W)(i - <V) (3) 

and: 

Vaa',PP' =U(1- S a/ 3)S aa ,Spf3> (4) 



We use V 1 and V 2 in the context of the GWA: the 
newly applied V 1 describes spin-flip scattering of an elec- 
tron upon interaction with an opposite spin electron, as 
shown schematically in Fig. HJa); V 2 does not allow 
such spin-flip scattering [see Fig. [TJb)] but has been 
found superior— ~ 27 ' 29 over the spin-independent version 
F 3 n , pp, = USaa'Spp' (V 3 is the form most commonly 
used in ab initio GW calculations of real materials). In 
the next section we show that V 1 outperforms V 2 when 
describing the low-energy properties of the SIAM. 

B. oGaW formalism 

We consider the SIAM in both equilibrium and 
nonequilibrium situations using a nonequilibrium Green's 
functions (NEGF) approach. Electron correlation ef- 
fects are obtained by solving self-consistently for the var- 
ious (retarded, lesser, etc ..) Green's functions and self- 
energies of the impurity. Details of our NEGF approach 
can be found in Refill; for simplicity we show in this sec- 
tion only time-ordered quantities relevant for the equilib- 
rium, zero temperature case. 

At half-filling (one electron on the impurity) the im- 
purity orbital energy is td — — (the zero of energy is 
chosen such that the chemical potential is fi = 0). Then, 
the impurity Green's function reads: 



ay ' u- [-U/2 + VH + Y, a (u)] +iA sgn(w) w 

Here, is the spin-dependent Hartree potential; one 
has Vj 1 = in the spin-flip case and V^ 1 = f/(n_ CT ) in 
the non-spin-flip case. 

The impurity self-energy £ is evaluated using the self- 
consistent GWA within the aGaW formalism 2 ^, which is 
based on the recent generalization of Hcdin's equations 
for quantum many-body systems with spin-dependent in- 
teractions. Following Ref3&, we decompose V in a basis 
formed by the Pauli matrices: 

Vaa',/30' = via'VlJVpp' ( 6 ) 

where repeated indices are summed and subscripts de- 
noted by capital letters (e.g. /, J) are indices represent- 
ing the four Pauli matrices a°,a x ,a v and a z . The V/j 
components associated with V 1 and V 2 are easily ob- 
tained from: 

VLf,fiP> =~2 (<a'°h' + °l a -°lp) ( 7 ) 

and 

Vaa',0P> = "2 ( a aa' a %' ~ a aa' a PP') ( 8 ) 

which shows that both V 1 and V 2 have contribution from 
the spin-channel. (In contrast, the spin- independent ver- 
sion V 3 has contribution only from the charge channel: 
^ 3 o = U, Vf^j^ = 0.) 
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Within the aGaW formalism the spin-dependent 
Hedin's equations read^S 

SaH = i J ^e* 0+E a^G^E^Wj^u E) 



W>j(w) = Vu + V ik P K l(u)W L j(uj) (9) 

where the irreducible polarization is evaluated within 
random-phase approximation (RPA), 



P u (lo) = -i 



J^^G p {E)aj a G a {E-uj) (10) 



The spin-flip and non-spin-flip GW flavors are obtained 
using definitions (J7J and (J5J for Vij in Hedin's Eq. (|9|). 
Using the following notation for the polarization bubble 
in spin space: 

/dE 
— G a {E)G p {E~Lu) (11) 
Z7T 

the final expression for S within the spin-flip GW flavor 
reads 

= -ij ge i0+ * G- a (E) 



III. RESULTS 

A well known characteristic of the SIAM is that its 
low-energy physics shows universal scaling behavior in 
the Kondo regime, i.e. near half-filling and for effective 
interaction strength: 



U 

7rA 



> 2 



(16) 



(a strong-interaction asymptotic regime is approached 
exponentially fast for u > 2). A signature of this aspect is 
that functions describing spectroscopic, thermodynamic, 
and transport quantities at low energies become inde- 
pendent of u when their energy argument is scaled by 
the Kondo temperature Tk- 



A. Scaling behavior of the inpurity density of 
states 

Figure [5] shows results at T = for the impurity den- 
sity of states: 



A(E) 



1 



\3tG(E)\ 



(17) 



U + U 2 P- a ^{u-E) 



1 + UP- a . a + VPa.-a + U 2 P- a . a P a ,- a 



(12) 



(above, we omitted for simplicity the energy variable on 
P showing in the denominator) . 

In the spin-flip case the expression for the self-energy 
of an electron with spin a involves a Green's function for 
an electron with opposite spin — cr, and similarly the po- 
larization bubble involves Green's functions of electrons 
with opposite spin. The same does not hold in the non- 
spin-flip case, where the self-energy expression reads^ 



^ e M+E GAE) 
Zn 



U 2 P„ 



,(w - E) 



1 - U 2 P^(uj - £)P_ CT> _ CT (a; - E) 



(13) 



In the absence of a magnetic field, the exact solutions 
for G and £ are non-magnetic, independent of spin. Sim- 
ilarly, we consider in this work only non-magnetic GW 
solutions in which case the expressions for E 1 and E 2 
simplify and read (we omit from now on spin indices): 



E i( w ) = -i\ d A e M+E G(E) £ 

' 2tt '\ + UP{uj-E) 



£ 2 M 



dE 
2tt 



JO+E 



G(E) 



U 2 P(uj - E) 
1-U 2 P 2 (uj - E) 



(14) 



(15) 



as function of energy divided by the half-width at half- 
maximum of the central peak (denoted -Ek)- Figure^a) 
shows results obtained within the spin- flip GW flavor, 
while Fig. [5{b) treats the non-spin-flip case. We see that 
GW captures the aformcntioned universality aspect: for 
u > 2, A(E) for various u are almost undistinguishable. 

From the numerically exact solution it is known that 
the low-energy part of the impurity density of states can 
be fitted very well with the following for m 30 ' 31 : 



A DS (E) 



1 

7rA 



'If 



i.A 



DS 



E + iA DS 



(18) 



with Ads/^k a constant determined by a via 
A» S (E K ) = 1/(2ttA). 

The heuristic form shown in Eq. (|18p is inspired by 
the Doniach-Sunjic law- 3 ^ that describes photoemission 
of core electrons by x-ray in metals. The role of the core 
electron in metals must be reflected in the SIAM by some 
sort of spin polarization cloud^ 3 . developed from electrons 
in the leads screening (at low T and low V) the unpaired 
impurity spin. 

For the numerically exact solution of the SIAM the 
best Doniach-Sunjic fit is realized for: 



CKexact — 0.5 



(19) 



This is believed to arise from a Fermi-edge singularity 34 
in the presence of a tt/2 phase-shift at the Fermi level. 
The scaling function associated with A(E) corresponding 
to the numerically exact solution is shown in Fig. [He) 
as the black dashed curve. 
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Also shown in Fig. [2{c) are the GW scaling func- 
tions representative of A(E) in the asymptotic regime: 
the thick solid (in red online) curve represents the spin- 
flip case and the thin solid (in green online) curve 
the non-spin-flip case. The Doniach-Sunjic fit works 
very well in the GW case as well (with an error 
ttA \A ds (E) - A(E) \ < 1), as demonstrated in Fig. HJc) 
by the black discontinuous curves. More precisely 



GW also displays universal scaling behavior for cr(T). 
The curves in red and green in Fig. [3] are representative 
of the spin-flip and non-spin-flip flavors respectively. We 
find that the Goldhaber-Gordon fit works for GW as well, 
as demonstrated by the black discontinuous curves in Fig. 
OH We also find that the value of the parameter s (sqw ~ 
0.11 in the spin-flip case, and sqw ~ 0-08 in the non-spin- 
flip case) is simply related to the a parameter through: 



0.26, spin-flip case 
0.19, non-spin-flip case. 



(20) 



We see that spin-flip GW does better than non-spin 
flip GW in describing the shape of the scaling function 
associated with A(E). However, the smallness of oqw 
relative to a exact indicates that both spin-dependent GW 
flavors strongly overestimate spin correlations between 
impurity and lead electrons. [In contrast, within spin- 
independent GW (or, more obviously, the Hartree-Fock 
approximation) which does not account for correlations 
in the spin channel, the Doniach-Sunjic fit is trivially 
realized for a — 1.1 



B. Scaling behavior of the linear-response 
conductivity 

Figure [3] shows the linear- response conductivity 



a(T) 



dl_ 

dV 



(21) 



v=o 



in units of the quantum conductance Go = 2e 2 /h and 
evaluated as: 



<y(T) 

Gq 



TTA 



dE^j^A(E) (22) 



where f(E) — l/[exp(E/T) + 1] is the usual statistical 
factor for a system of electrons in equilibrium at finite 
temperature T and /i = 0. 

Universal curves are obtained in the asymptotic regime 
w > 2 by plotting cr(T) as a function of T ' /Tk where Tk 
is the Kondo temperature. There are several definitions 
of Tk existing in the literature, given that the Kondo 
scale can be defined in the asymptotic regime only up 
to a multiplicative factor. Here, we define Tk as the 
half-width at half-maximum of the linear response con- 
ductivity: <j(Tk) = 0.5Go- 

The numerically exact, universal results for a(T) are 
represented in Fig. |3]by the black solid curve, obtained 
here from the empirical Goldhaber-Gordon forrn^: 



_GG 



(T) 



T 



K 



Go 



T 2 + T'£ 



(23) 



where T K = Tk / ^f2 1 l s — 1 and the parameter s takes the 
value Sexact = 0.22 when the Goldhaber-Gordon form is 
used to fit NRG results^. 



Sexact 



C^cxact 



(24) 



For T <C Tk, all curves in Fig. [3] show Fermi- liquid 
behavior, 1 - cr(T)/G ~ {T /Tk) 2 ■ Due to the smallness 
of s (or equivalcntly of a), a drops too fast with T/Tk 
within GW, even though the spin-flip flavor improves sig- 
nificantly over the non-spin-flip one. For T Tk, o 
drops too slowly, as expected since GW cannot capture 
the Coulomb blockade regime unless it artificially breaks 
spin symmetry. 



C. Scaling behavior of the differential conductivity 

Figure 2] shows G W results for the differential conduc- 
tivity at T = 0: 



T=0 



calculated from: 



— — = irAA — 
Go V2 



x ; y/2 JT ^dA{E) 
ir A I dE- 

-V/2 



dV 



(25) 



(26) 



and plotted for u > 2 as a function of a symmetrically- 
applied bias voltage V divided by a Kondo scale Vk de- 
fined as (t(Vk) = 0.5Go- 

The GW results for a(V) show again universal scaling 
behavior, i.e. there is just one representative curve for 
any u > 2. The curves in red (with circles) and green 
(with squares) shown in Fig. U correspond to the spin- 
flip and non-spin-flip flavors respectively. (We are not 
aware of numerically exact results that can be compared 
to our GW data in Fig. gj) 

We find that the GW calculated <r(V) can be fitted 
very well by the Goldhaber-Gordon form a GG (V), as 
demonstrated by the discontinuous black curves in Fig. 
21 The value of s is ~ 5% larger than the corresponding 
value in the case of cr(T), indicating that the effect of 
a low bias voltage is quite similar (but not exactly the 
same^i) to that of an effective temperature. 



D. Kondo scale 

We now turn our attention to the Kondo scale that 
joes inside the scaling functions and compare in Fig. [S] 
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the GW and analytic (asymptotically exact) Kondo tem- 
perature for a broad range of the interaction strength in 
the asymptotic regime. 

The Kondo temperature of the exact solution is well 
known 3 - to decrease exponentially with u, 



Tk 
A 



(27) 



and it is shown as the black solid line in Fig. [5j Also 
shown by the curves in red (with circles) and green (with 
squares) are the GW results obtained within the spin- 
flip and non-spin-flip flavor respectively. Importantly, we 
find that the GW Kondo scale depends algebraically on 



K 

A 



(28) 



as indicated by the black discontinuous curves in Fig. 
03 The value of the parameter /3 that controls the GW 
Kondo scale is: 



3.2, spin-flip case 
5.1, non-spin- flip case 



(29) 



The large value of j3 implies that the GW Kondo scale 
decreases fast with increasing interaction strength. For 
this reason, finding a non-magnetic GW solution for u > 
2 requires a very fine sampling of the energy axis near 
Ep, that we achieve using a discretized logarithmic scale 
~ ±A~™ (A as small as 1.005 and n as large as 4000 
were used for the largest u considered). In the Appendix 
we provide some simple arguments that relate /3 to the 
parameter a that controls the GW scaling functions. 

We also note from Fig. [5] that spin-flip GW performs 
better than non-spin-flip GW for the most physically rel- 
evant u- values (i.e., it < 6). At strong enough interaction 
both GW flavors are guaranteed to severely overestimate 
the Kondo scale due to their algebraic u dependence. 

We conclude this section with the relationship between 
the Kondo scales Ek, Tk and Vk that go inside the 
scaling functions associated with A(E), <j(T) and a(V). 
Within both spin-dependent GW flavors considered we 
find: 



GW 



i.2 r, 



GW 



K 



and: 



V, 



GW 



K 



1.5 T, 



GW 



A" 



(30) 



(31) 



while for the numerically exact solution, NRG yields (see 
Ref42): 



2.3 Tf 



(32) 



E. Fermi-liquid properties and the linear 
coefficient of the specific heat 

The Fermi-liquid properties manifest themselves at 
very low energies: w, T, V <C Tk- Here, is character- 
ized by quadratic behavior, and in the asymptotic regime 
the impurity self-energy has the following low-order ex- 
pansion 



E(w,T, V) 



~Z 



iC 
2AZ 2 



2 rp'2 



it T 



(33) 



For the exact solution, the form shown in Eq. (|33|) can be 
derived via the Ward identity, and one find o 12 ' 41 " — C = 
1. Within the GWA the Ward identity is not satisfied 44 
and we find instead C ~ 2.6 in the spin- flip case and 
C ~ 3.5 in the non-spin-flip case. 

The Fermi-liquid behavior is dictated by the inverse 
renormalization factor 



Z~ 



1 



as 

duj 



(34) 



uj,T,V=0 



which is directly related to the T-linear coefficient of the 
impurity heat capacity 



1 dE imp 
7 = lim — - 

' t^o T dT 



(35) 



where i?i m p is the average energy associated with the 
presence of the impurit} ^ 27 ' 45 . Indeed, normalizing 7 to 
its value in the absence of Coulomb interactions 



3A 
2^ 



7 



(36) 



one obtains that 7 is enhanced by Coulomb interactions 
by a factor equal to Z~ Y [we have checked this statement 
in the GW case by directly evaluating the right-hand side 
ofEq. {35])], 



(37) 



Figure [6] shows the ratio 7Gw/7exact as a function of 
u, for both GW flavors considered in this work ( 7 cxa ct is 
evaluated according to Ref4£). One can see that within 
non-spin-flip GW, 7 is largely overestimated over the en- 
tire u-range considered; the spin-flip GW version signifi- 
cantly improves 7 for u < 9, but the description remains 
unsatisfactory. In the limit of very large u both GW 
flavors are bound to underestimate 7 since universality 
implies Z~ x - A/T K . 

To analyze more in detail the reasons behind this in- 
correct description, we note that within spin-dependent 
GW 



J GW 



dE A(E)—mV(E) 
dE 



A(0)W(0), (38) 



(We are not aware of numerically exact results for 
Vk/T k .) 



where we used the facts that within RPA (i) the energy 
dependence of W is dictated by a pole situated on the 
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imaginary axi; 



,47 



W(w -> 0) 



IUJ — (jJsf 



(39) 



at the so-called spin- fluctuation frequency losf that is 
directly related to the static polarization bubble P(0)r^ 



7rA 2 



1 + UP(0) 



and (ii) we find that wsf <C -Ek- 
To this end one can write: 



Z, 



GW 



A 

WSF 



(40) 



(41) 



and conclude that the incorrect description of Zq W is 
intimately related to the incorrect description of the spin- 
fluctuation frequency u>sf [or of the static polarization 
bubble P(0)] within RPA. 

The above analysis suggests that vertex corrections be- 
yond RPA in the irreducible polarization P are important 
for a correct description of spin fluctuation effects and of 
the Kondo peak. This complements a finding of previ- 
ous work2£, namely that vertex corrections beyond GW 
in the electron self-energy E are critical for a correct de- 
scription of charge fluctuation effects and of the Hubbard 
bands. 



IV. SUMMARY AND CONCLUSIONS 

We have found within GW that the low-energy prop- 
erties of the SIAM display universal scaling behavior in 
the Kondo regime. This is remarkable if one thinks of 
GW as a theory based on perturbation expansion about 
the impurity Coulomb interaction U, but perhaps not 
very surprising given that GW corresponds to a partial 
resummation of an infinite set of such diagrams. 

We have shown that the scaling functions describing 
quantities such as the impurity density of states or the 
linear-response conductivity have a similar form within 
GW as in the numerically exact solution, but they are 
characterized by a different parameter value. The de- 
viation of this value from the correct one indicates the 
extent at which spin-dependent GW overestimates spin 
correlations between the impurity spin and the electrons 
in the leads. We have also found that the Kondo scale 
that goes inside the GW scaling functions depends al- 
gebraically on the interaction strength, as opposed to 
the exponential dependence shown by the exact solution. 
Both the parameter value characterizing the scaling func- 
tions and the Kondo scale at physically relevant interac- 
tion strengths are best described within the spin-flip GW 
flavor. However, none of the GW flavors describes satis- 
factorily the Kondo regime. 

The importance of our analysis should be two-fold. 
First, it clearly quantifies the extent at which various GW 
flavors fail to describe satisfactorily the Kondo regime 



of the Anderson model. Second, it suggests that a sys- 
tematic study of the low-energy properties of the An- 
derson model is possible within a perturbative approach 
based on expansion about the screened Coulomb interac- 
tion W: with increasing order of perturbation expansion 
about W, the parameters controlling the scaling func- 
tions should converge toward the exact values, while the 
Kondo scale should improve overall with best accuracy 
expected at lower u. Along this line, extending Hedin's 
approach up to second order expansion about W ap- 
pears naturally as a path towards improvement beyond 
GW/RPA. 
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Appendix: Simple arguments relating the 
parameters that control the GW scaling functions 
and the Kondo scale. 



Within GW, we find that the low- and high-energy 
parts of the impurity (retarded) Green's function have 
the following forms in the strong interaction regime: 



G(E) 




(A-l) 



where E x > T K and E 2 <.U. 

The low-energy (\E\ < E%) behavior of G(E) is consis- 
tent with the Doniach-Sunjic form for the impurity spec- 
tral function [see Eq. (fT5)) ]. With good approximation 
we find that E\ can be pushed up to energies of order A. 

The high-energy (\E\ > E 2 ) behavior of G(E) is con- 
sistent with (if not a consequence of) the fact that within 
RPA one has (to leading order in l/U) 



P(0) = - 



1 



{G 2 (P)}«- 



(A-2) 



a result intimately related to the smallness of the GW 
Kondo scale [see Eq. (|4"0")) ]. With relative good approxi- 
mation, Ei can be pushed down to energies of order A. 

Let us assume for the moment that the low- and high- 
energy behaviors shown in Eq. (IA-1|) are exactly valid 
in the intermediate-energy regime near E ss A ^> Arjs- 
Then, consistency in the P-dependence of G at fixed U 
implies f(x) ~ x~ a for x <C 1, while consistency in the 
[/-dependence of G at fixed E yields 



1 



(A-3) 
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Eq. (|A-ip hold only approximately near E ~ A and 
the relation between the parameter (3 controling the GW 
Kondo scale and the parameter a controling the GW 
scaling functions is in practice 

one obtains within the above assumption that: 

P^l/a-1. (A-5) P>l/a-l. (A-6) 

In reality, the low- and high-energy forms shown in 



Since [see Eqs. (18]) and (l28l) ] 

Ads T k 



A A ^ 
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FIG. 1: Diagrammatic representation in spin space of the 
local two-particle Coulomb interaction V: (a) spin-flip case, 
v aa',pp' = -U{l-8 a p){l-5 aa >)(l-S ppl ), and (b) non-spin- 
flip case, V a 2 a , i/3(3 , = 17(1 - 5 a p)8 aa ,5pp> . 
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FIG. 2: The impurity density of states at low energy and 
T = 0, V = 0. (a) Results within the spin-flip flavor of 
spin-dependent GW, for several values of the effective inter- 
action strength u. (b) Results within the non-spin-flip flavor 
of spin-dependent GW, for the same u values as in (a), (c) 
Comparison between spin-dependent GW and the heuristic 
Doniach-Sunjic (DS) fit to NRG results in the asymptotic 
regime. 



10 



1 

0.9 
0.8 
0.7 



O 0.6 
to 

c 0.5 

e °- 4 

e 0.3 
0.2 
0.1 







-O-GW spin-flip, u>2 




-D— GW non-spin-flip, u>2 




GG fit to NRG, s=0.22 - 




GG fit, s=0.11 


*** ^ \ 


GG fit, s=0.08 


: X 


x. ** 







0.01 



0.1 1 

T/T, 



10 



100 



K 



FIG. 3: Linear-response conductivity in the asymptotic 
regime at low temperature, calculated within spin-dependent 
GW and from the empirical Goldhaber-Gordon (GG) fit to 
NRG results. 
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FIG. 4: Differential conductivity in the asymptotic regime at 
low bias voltage and T = 0, calculated within spin-dependent 
GW. 
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FIG. 5: Comparison between the Kondo temperature cal- 
culated within spin-dependent GW and the one from the 
(asymptotically exact) analytic solution (Ref.— ) for a large 
range of the effective interaction strength u in the asymptotic 
regime. 




FIG. 6: Ratio of the impurity heat capacity calculated within 
spin-dependent GW to the one calculated exactly (by the 
Bethe-Ansatz method) for a large range of the effective in- 
teraction strength u. 



